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SUMMARY 


Turbulent  and  precipitation  conditions  in  the  real  atmosphere  are 
considered  in  their  effect  on  the  line-of-sight  propagation  of  electromagnetic 
waves.  It  is  found  that  tropospheric  scattering  and  phase  scintillation 
mechanisms  place  an  interesting  limit  on  the  accuracy  with  which  one  may 
determine  the  angular  position  of  a  distant  radiator.  Superposition  of  the 
same  perturbed  signals  gives  rise  to  waveform  distortion  in  communication 
systems.  Theoretical  expressions  for  angular  power  distributions  and  phase 
error  spectra  are  derived  from  phenomenological  atmospheric  models.  Ray 
bending  by  cloud  formations  and  rainstorms  is  also  computed  for  representative 
weather  conditions. 
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I.  INTRODUCTION 


The  present  report  summarizes  the  results  of  an  investigat  on  initiated 
to  answer  *be  simple  question:  "With  what  precision  can  one  propagate  micro- 
wave  electromagnetic  energy  along  a  Tine -of-sipht  path  in  the  real  atmosphere*^" 
Such  a  question  arises  naturally  if  one  asks  for  the  ult  mate  accuracy  of  radio 
location  measurements  or  ultra -high  speed  communication  techniques.  One 
car  imagine  a  variety  of  scattering  and  phase  shifting  mechanisms  which 
might  operate  in  the  atmosphere  to  produce  angle  of  arrival  problems,  multi - 
path  propagation  modes,  signal  phase  scintillation,  ar.d  broadening  of 
received  energy  patterns. 

This  study  was  fortunate  in  that  considerable  attention  has  recently 
been  focused  on  the  atmosphere  in  connection  with  beyond-linc-of-sight 
(scatter)  communication  systems.  A  theory  of  turbulence  scattering  has  been 
proposed^ ^  ^ to  account  for  observed  propagation  anomolies.  There  has 
subseq  itrtly  been  initiated  a  thoroughgoing  examination  of  the  atmosphere 
and  its  variations.  In  addition  to  the  usual  meteorological  disturbances 
(snow,  rain,  clouds,  and  dust),  an  atmospheric  fine  structure  extending  from 
the  very  surface  of  the  earth  to  ionospheric  heights  has  thereby  been  uncovered. 
This  structure  is  measured  as  time  and  space  varying  concentrations  of  index 
of  refraction  in  the  troposphere  and  free  electron  densities  in  the  ionosphere. 


E.  Gordon,  Proc.  I.R.E. 
I  R.E.  ,  35,  453  (1947K 


41  ,  284  (1950). 
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The  present  study  h^is  exploited  phenomenological  models  for  turbulence 
and  precipitation  scattering  which  have  been  fitted  to  such  experiments.  Effects 
of  clouds  ard  elevated  layers  are  computed  from  estimated  departures  of  the 
index  of  refraction  from  its  median  values.  In  general,  no  attempt  has  been 
made  to  determine  the  effects  of  signal  attenuation  *  ducting  or  uniform 
horizontal  layers.  The  essential  output  of  the  following  calculation  is*  (1)  RMS 
phase  error  and  scintillation  spectra  for  single  rays  and  for  linear  superposi¬ 
tion  of  multipath  signals,  and  (2)  angle  of  arrival  variations  for  single  rays 
and  angular  distributions  for  the  scattered  energy.  No  account  is  taken  in 
this  report  of  the  stratified  nature  of  the  atmosphere,  only  deviations  from 
the  mean  are  considered.  These  results  are  computed  for  one  receiving 
antenna  only,  so  that  in  order  to  predict  performance  limits,  one  must  employ 
the  characteristics  of  the  particilar  system  contemplated.  This  study 
indicates  two  points  rather  clearly:  (1)  the  atmosphere  can  be  a  limiting 
factor  in  very  accurate  radio  location  schemes,  and  (2)  at.tendcnt  phase 
variations  imply  a  real  distortion  limit  for  ultra-high  speed  communication 
systems . 

II.  TROPOSPHERIC  TURBULENCE 
The  troposphere  is  in  continuous  turbulent  motion,  and  is  therefore 
neither  uniformly  stratified  nor  homogeneously  mixed.  At  a  given  instant 
there  are  both  horizontal  and  vertical  fluctuations  in  the  temperature, 
pressure,  and  humidity  about  their  respective  means.  One  may  think  of 

Except  insofar  as  this  is  necessary  to  secure  convergence  (and  thus  proper 
physical  description)  for  wide-angle  acceptance  receiving  antennae. 


-  3  - 


these  micromctpor ological  disturbances  as  a  turbulent  boundary  layer 
phenomenon  associated  with  motion  of  the  winds  over  the  earth's  rough 
surface.  A  dela  ted  eddy  decay  scheme  has  been  suggested  which  imagines 
these  large  scale  wind -produced  turbulence  clusters  to  subdivide  until  their 
size  and  energy  are  eventually  absorbed  in  moler  iilar  friction.  We 

shall  base  our  development  on  a  purely  phenomenological  description  of  these 
tropospheric  variations.  The  statistical  properties  of  the  troposphere's 
dielectric  constant,  and  hence  its  index  of  refraction,  are  given  by  the 
normalized  space  correlation  function 


p{r,2)  = 


Ae(r) 


(2.  1) 


whore  r  vector  distance  between  points  (1)  and  (2).  Ae  is  the 

var  lation  of  the  dielectric  constant  about  its  mean.  This  function  defines 
the  (normalized)  cross  correlation  between  simultaneous  records  of 
dielectric  constant  taken  at  two  points  a  distance  r  apart.  The  atmosphere 
is  assumea  to  be  isotropic  and  homogeneous,  so  that  c(r^  )  is  a  function 
neither  of  the  direction  nor  the  end  points  of  . 

Limited  tests^^^’^^^  tend  to  support  the  exponential  correlation 


^  F .  Villars  and  V .  F.  Weisskopf,  Phys.  Rev.,  94,  232  (1954).  They  have 
given  a  fuller  theory  of  turbulence  scattering  b«asTB  on  the  Navier -Stokes 
equation.  Their  paper  also  contains  rcfeVenccs  to  the  previous  work  of 
Heisenberg,  Batchelor,  and  others. 

MacCreddy,  Jr.,  Jour,  of  Meteor.  ,  10,  434  (1953). 

(5) 

J.R.  Gerhardt,  C.M.  Crain,  and  H.  W .  Smith,  Jour,  of  Meteorology,  9, 
299  (1952). 

Birnbaum,  Phys.  Rev.,  82,  110  (1951). 
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P(ri2)  =  e 


(2.2) 


where  is  a  measure  of  the  scale  of  the  turbulence, 
o 

Detailed  measurements  have  shown  that  the  isotropic  assumption  is 

not  strictly  true;  the  scale  length-/^  differs  in  horizontal  and  vertical 

(71 

directions,  at  least  near  the  ground.  Direct  measurement'  '  has  also  shown 

that  this  scale  increases  somewhat  with  altitude.  Var  iations  in  the  measure > 

(71 

ments  of  scale  length  with  location  and  time  of  day  ai  ?  so  large  that  it  was 
not  considered  worthwhile  to  use  anything  more  sophisticated  than  Equation 

(2.2)W. 

Using  Eq.  (2.2),  Booker  and  Gordon^^^  have  shown  that  the  power 
scattered  per  unit  solid  angle,  per  unit  incident  power  density,  and  per  unit 
volume  is  given  by 


0(9, 


(2.3) 


\  denotes  the  wavelength  of  the  incident  radiation  and  tL  is  the  angle  between 
the  incident  electric  field  vector  and  the  direction  of  scattering.  6  is  the 
angle  included  between  the  direction  of  incidence  and  scattering  as  shown  in 
Fig.  1. 


T71 - 

'  'C.M.  Crain,  A.  W. Straiton,  andC.E.  VonRosenberg,  Trans.  I.R.E.,  Vol. 
AP-1,  No.  2,  October,  1953. 

L.  Staras ,  Jour.  Appl.  Phys.  ,  23,  1152  (1952),  has  investigated  the 
effect  of  other  assumptions  for  p(r)  on  scattering  phenomena. 
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Figure  1 

Coordinates  for  Turbulence  Scattering  of 
Vertically  Incident  Radiation 
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Figure  2 

Multipath  Propagation  by  Scattering 


A.  Scattered  Wave 


We  shall  consider  first  the  energy  scattered  into  a  single  antenna  by 
off-axis  turbulence  clusters.  The  situation  is  pictured  in  Fig.  (2)  where  the 
multipath  propagation  mode  is  clearly  exhibited.  The  scattered  waves  give 
rise  to  an  angular  distribution  of  received  power  about  the  antenna's  pointing 
axis,  as  indicated.  The  ratio  of  the  power  received  by  the  antenna  along  the 
multipath  (R  +  r)  to  that  incident  along  the  line-of-sight  (Rq)  ^or  a  scattering 
volume  dv  is 


dP  R^cr(e.») 

— i  =  - dv 


(2.4) 


o 

where  a  denotes  the  scattering  cross  section  per  unit  volume  per  unit  solid 
angle.  One  may  set  =  R  for  reasonably  (  <[4°)  narrow  beams  and  long 
range  (>1000m)  transmission. 


dP 


2 


(2.4a) 


The  p»ath  difference  between  the  direct  and  multipath  ray  is  simply 


X  =  R  +  r  -  R^  .  (2.5) 

If  one  integrates  (2.4a)  over  all  volume  elements  satisfying  this  condition, 
the  result  may  be  interpreted  as  a  "differential  path  difference  power 
density,"  dP  =  Q(x)dx.  The  ratio  of  the  scattered-to-direct  power  received 
with  path  difference  x  is  thus: 


Q(x) 


<t(0.  ^) 

2 

r 


6{x  +  R^ 


R 


r)  , 


(2.6) 
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P 
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where  6(z)  ig  the  familiar  Dirac  delta  function. 

(1)  Phase  Delay 

Equation  (2.6)  may  now  be  applied  to  predict  an  RMS  phase  error  in  a 
sinusoidal  voltage  at  the  receiving  antenna.  The  anomolous  propagation 
produces  a  linear  superposition  of  scattered  (and  thus  delayed)  sinusoidal 
signals,  as  pictured  in  Fig.  (3a).  The  resultant  phase  error  (a)  so  produced 
maybe  inferred  from  the  vector  voltage  diagram  of  Fig.  (3b).  The  lag  angle 
between  the  primary  voltage  (assumed  for  the  moment  to  be  of  constant 
phase)  and  scattered  voltage  E(x)  is  0  =  2itx/X.  .  The  desired  phase  error  a 
is  computed  as: 


(2.7) 


in  the  approximation,  E(x)  <^E^.  Introducing  Eq.  (2.  6),  the  mean  squared 
error  becomes: 


/ r  .  2/2wx\Pj  O'(0.^)cr  .  D  D 

<  Q  >  =  \  dx  sin  { - )  \  dv  — i 6  X  R  -  R  - 

^  /  Jo  VJJ  L  ° 


(2.8) 


The  primed  volume  integration  restricts  that  process  to  the  admittance  angle 
(i.e.,  beamwidtl)  of  the  receiving  antenna.  This  step  is  a  crucial  one,  since 
widening  that  angle  indefinitely  would  only  serve  to  include  all  scattering 
elements  within  sight,  and  so  produce  a  d.vergent  answer.  The  spherical 
shape  of  the  earth  and  limited  (height)  extent  of  the  turbulence  reduces  such 
an  accumulation  in  actual  experiments.  One  must  also  consider  the 
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Figure  3 

Superposition  of  Multipath  Signals 
(a)  Principal  Wave  Plus  Delayed-Scattered  Signal 
(b)  Vector  Voltage  Diagram 


attenuation  of  the  scattered  waves  for  a  complete  treatment.  We  shall  avoid 
such  cOiisiderations  by  exploiting  the  characteristically  small  beamwidths 
of  radar  rece.vers  as  a  natural  cut-off. 


To  compute  tie  integrals  appearing  in  Eq.(2  8),  we  erect  a  spherical 
coordinate  system  on  tie  receiving  antenna's  po  ntlng  direction,  as  shown  in 
Fig.  2.  If  p  is  the  antenna  beamwidth  and  L  the  transmission  range  through 
the  troposphere, 

2/2.,x\r^_ 


^  dx  Bin^f  dr  ^  d0  einS  ^  df  (T(6,0)6rx  -  r(l  -  co80)^ 

^  '  o  \  X  /’"o  o  o  L  J 


(2.9) 


We  must  now  express  <j{Q  d'  in  terms  of  0  and  p  . 

2Tr./  ^ 


d(e,»)  = 


(  °\  r. 

^  2\  -12  \  /  ^ 

CAN  >  10  \  X  ^  L 


^(6)  co8^(0  -  y^) 


,2ttC\2 

1  +  2  \^—^J  ( ‘  * 


(2.10) 


The  definition  of  the  dielectric  constant  in  terms  of  the  index  of  refraction 
(i.e.,  e  =  n^)  and  the  N-unit  notation,  AN  =  (n  -  1)10^  are  here  utilized.  All 
but  one  of  the  integrations  in  Eq.(2.9)  may  be  performed  analytically  if  one 
disregards  altitude  variations  in/^  and  AN^  —  a  procedure  which  is 
certainly  consistent  with  our  other  approximations. 


^  -  <AN^^ 


1  - 


sin(u)'' 


u 


\  L 


where  .  . 

U  =  mil  (1  -  cosP)  . 
./  o 


(2.11) 
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oo‘s 


The  ratio 


is  plotted  in  Fig.  (4)  as  a  function  of  2tt^  /\ 


for  various 


L//  . 
o 

(2)  Angular  D’str.bution 

The  angular  distribution  of  scattered  power  about  the  receiving 
antenna's  pointing  direction  may  be  extracted  from  Eq.  (2.  6)  by  integrating 
over  all  path  differences  x.  Again  using  an  antenna-centered  spherical 
coordinate  system. 


cosp 


L 

cos0)j  d^<T(0,^)  , 

o  o 


(2.  12) 


where  or (9, 0)  is  given  by  Eq.  (2.10).  The  angular  power  di stribution  is 
recognized  as  the  d(cos0)  integrand  in  Eq.(2.12),  so  that 


1  dP 
s 

d{cos0) 


2L 

T* 

o 


<AN^ 


10 


,-12 


1  +  cos  (0) 


1  +  2 


/2tT/\ 


(1  -  COS0) 


(2.  13) 


where  constancy  of  the  scattering  parameter s ,  and  <s^N  ,  over  the 

transmission  path  has  again  been  assumed.  The  distribution,  Eq.  (2.  13),  is 
plotted  in  Fig.  (5)  as  a  function  of  0  for  various  2tt.^/\  .  In  the  region 

0  ^6^ — - —  2:  denominator  is  sensibly  unity  and  the  density  equals 

2tt,^ 

o 


4 


Zn  C 
o 


For  larger  angles,  it  falls  off  rapidly,  reaching  1/4  at  right 


angles  to  the  beam. 

The  result,  Eq.  (2.  13),  may  be  further  re  fined  by  including  the  effect 

of  attenuation  of  the  scattered  signals  in  reaching  the  antenna  from  the 

(9) 

various  volume  elements.  LaGrone,  Benson,  and  Straiton  have  computed 

the  energy  which  is  removed  from  a  beam  per  mile  with  the  same  theory  of 

tropospheric  turbulence  scattering  used  here.  Their  result  is 

'2Tri  \ 


db 


(2.  14) 


mile  V  X  /  '  ^  ^  / 

In  evaluating  Eq.(2.  12)  one  should  thus  include  a  range -dependent  absorption 
factor 


1  L  27t 

^  C  d(cos0)  C  dr  C  d0 
^cosp  - 


ct(0,^)  e 


o  ~o 

_  r|700  —  <AN^  10" 
X 


12 


X  . 


(2.  15) 


The  attenuation  effect  is  almost  completely  negligible  in  this  case,  since 
forward  scattering  is  so  favored  by  the  cross  section,  Eq.  (2.  3). 

B.  The  Direct  Wave 

The  preceding  discussion  was  concerned  only  with  the  scattered  wave. 
As  in  any  scattering  problem,  the  resulting  field  at  a  point  in  space  is  com¬ 
posed  of  both  the  direct  and  tcattered  wave.  In  the  present  instance,  however, 
the  direct  wave  cannot  be  written  as  a  simple  sinusoidal  plane  wave  since  it 
has  traversed  a  region  of  random  and  variable  delay. 


Jour.  Appl.  Physics,  22,  672  (1951), 
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( 1 )  Phase  Delay 


To  compute  the  phase  delay  suffered  in  traversing  a  length  L  of 

turbulence,  we  shall  write  the  received  electric  field  (Ej^)  in  exponential 

(10) 


form. 


=  E.J,  exp 


iZvf  ^ 
o 


d  I 

V 


(2.  16) 


E.p  is  the  transmitted  field  strength  and  f  the  frequency.  The  velocity  of 
propagation  (v^)  depends  on  the  dielectric  constant  and  so  is  a  function  of 
position  along  the  path.  If  we  separate  a  random  portion  from  the  dielectric 


constant,  e  =  +  Ae  ,  we  can  rewrite  Eq.  (2.  16)  as. 


Ej^  =  E.J,  exp  i 


inL  +  i  r  d^ 

X.  X  ^o 


Ae(Je)‘ 


o  -> 


(2.  17) 


The  random  part  of  the  phase  angle  is  represented  by  the  exponent's  second 
term.  Calling  this  a  and  taking  the  mean  square, 

^2  L  L 

<a^  =  d/jJ"  d/^  <  - 1- - ^  >  .  (2.18) 


The  integrand  is  just  ^  P(^j  ~  ^2}  definition  Eq.  (2.  1),  so  that 


one  has 


<•'>  ■ 


o  o  L  J 

(continued  on  following  page) 


See  also  J.  Feinstein,  Trans.  I.R.E.,  PGAP,  Vol.  AP-2,  63  (1954). 


L 


(2.  !9) 


for  the  mean  square  phase  deviation  of  the  direct  wave. 

To  determine  the  time  spectrum  for  this  turbulence,  consider  a  mass 
of  turbulent  air  mov  ng  with  an  average  velocity  v  Take  a  set  of  coordinates 
also  moving  with  the  mass  at  speed  v).  Measurements  of  fluctuations 

in  this  moving  system  will  give  data  from  which  the  time  spectrum  can  be 
calculated.  Unfortunately,  such  data  are  not  available  and  in  their  absence 
an  approximate  treatment  must  be  used.  In  the  following  we  shall  calculate 
the  time  spectrum  produced  by  a  line -of- sight,  rota *ing  through  a  turbulent 
air  mass  containing  variations  which  are  not  changing  in  time  relative  to  one 
another.  This  is  equivalent  ta  assuming  that  the  line -of  -  sight  moves  through 
a  characteristic  length  in  a  time  short  compared  to  the  time  of  significant 
change  in  the  turbulence  pattern.  In  Part  IV  we  shall  consider  the  analogous 
problem  of  the  line -of-sight  moving  parallel  to  itself.  That  this  procedure 
will  yield  useful  results  is  suggested  by  some  experiments  performed  by 
Straiton  and  Smith. 

A  rotating  line-of-sight  sweeps  through  an  angle  6  =  t0  in  time  t, 
where  0  is  the  angular  rate.  (See  Fig. 6.)  Defining  a  as  before,  one  can  write 


<;^Ae*(x^,  yj)  Afe(x^, 


(2.  20) 


W. 


Straiton  and  H. 


W.  Smith,  Proc.  I.R.E.,  38,  825  (1950). 
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Figure  6 

Geometry  for  Rotating  Line -of -Sight 


Figur '  7 

Power  Spectrum  of  Phase  for  Rotating  Line -of-Sight 


where  the  first  integral  is  taken  along  path  at  time  t,  and  the  second  along 
line  ^2  t  seconds  later.  Using  Eq.(2.2)  and  measuring  distance  along  the 
line -of-sigbt  from  the  origin,  £q.(2.20)  can  be  written 

'  '  V  ^  o  o 


exp 


i: 


^1  +^2  ■ 


(2.21) 


It  is  convenient  to  assume  that  ?  L2  •  (i.e.,  that  the  path  is  long  compared 
to  any  charges  in  length),  so  that  Eq.(2.21)  can  be  changed  to  polar  coordinates 
and  integrated  once  to  yield: 


^  ^  lO'^^i^  it/4  f 

<^a(t)a{t  +  t)^  =  - ^ ^ d^<J- 


exp 


L^l  -  sin20  cos6 
cos  0 


-  ain2^  cos6 

L^l  -  sin2^  cos6 


exp 


/ 


o 


cos^  -  8in2^  COS0 


1  -  sin2^  COS0 


(2.  22) 


The  first  term  is  readily  integrated  again.  If  L/^^»l,  the  last  two  terms 
make  significant  contributions  only  for  0«tt/2  and  ^  =  Tr/4  .  One  may  put 


sin2^ 


^  J  N 


_  TT  \  2 


cos^  ^—L 
-\T 


(2  23) 


and  with  the  substitution 


0  = 


1  -  CO80  j  2  1  .  TT 

-  ly  -  1  +  -  , 

2  CO80  ^  4 


(2  24) 
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the  second  term  in  Eq.(2.22)  becomes 


tt/4 

3  dp 


.  I  L-\l  -  sin2o  COS0 
L  exp  I - - 

COBp 

o  ^ 
co80-^l  -  sin2^  cosG. 


/^cosS 

o 


/d. 
1 


exp 


2( i  -  cos6)  Ly 


/ 


4 


y  -  1 


(2.25) 


where  _ 

I  ^  A  ' 

Y  =  |_!L_£02«_  f  1  . 

8(1  -  cosO) 

When  1  and  0  <C<jt/Z,  one  can  put  Y 

thus  becomes 


(2.26) 


=  00  with  small  error.  The  integral 


exp 


T  ® 

-r  ■ 


dy 


J  L  /IGIL' 
- K 


(2.27) 


where  K  is  a  modified  Bessel  function  of  the  second  kind.  The  last  term  can 
o 

be  treated  in  a  similar  fashion. 


1  r® 

FT  j 

e  *^1 


exp 


J 


yjy  -  1 


1  ^00  oo  expR-xiy)] 

"3^i=r' 


2  9 


ttL 


•^r\ 


-X’r 


(2.  28) 


L  denotes  the  modified  Struve  function.  Combining  all  of  these  results,  one 
finds : 
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where  0  =  0T  .  When  t  =  0,  Eq.  (2. 29)  reduces  to  Eq.  (2.  1 9). 


The  spectrum  may  be  found  by  taking  the  cosine  transform  of 
Eq.  (2.29) 


W(f)  =  \  dr  cos2TrfT  <^Q(t)Q(t  +  t)^  . 


(2.  30) 


-00 


upon  substituting  Eq.  (2.  29)  into  Eq.  (2.  30),  the  second  term  may  be  integrated 
directly.  The  first  term  may  be  simplified  by  using  the  identity 

00 


z 


-  2  ™ 

1--  "  K^(x)dx  =  K^(z)L_j(2)  +  K_j(z)L^(z) 

_  TT  J 


The  integral  to  be  evaluated  can  then  be  written 


oo  cos  OJ.U 

r  du  ‘ 

‘  2 

1  -  - 

r  K  (x)  dx 

"o  u 

TT 

U  J 

(2.31) 


where  cjj  =  2nf-^/ L0  .  Equation  (2.  31)  may  be  integrated  once  by  parts  to  give 


4/  tr 
o 

0L 


oo 

dll  K^(u) 
o 


U>  jU 


cos  t 

- dt  .  . 


t 


(2. 32) 


By  changing  the  variable  in  the  second  integral  to  y  =  t/wjU,  the  order 
of  integration  may  be  reversed  whereupon  the  remaining  integrations  may  be 
performed  directly.  The  complete  result  is 
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W(f)=:  - log 

2* 

X  6 


1  +  /I  +  w 


1  +  w. 


(2.33) 


This  is  plotted  on  Fig.  7. 

(2)  Fluctua  tions  in  Angle  of  Arrival 

In  addition  to  the  phase  delay  suffered  in  traversing  the  troposphere, 

the  direct  ray  will  exhibit  a  random  variation  in  angle  of  arrival  at  the  receiver. 

(121 

This  effect  is  analogous  to  the  angular  scintillations  of  star  images  '  and  is 

produced  by  j  radients  In  the  index  of  refraction  at  right  angles  to  the  ray  path. 

The  equation  governing  propagation  of  a  ray  through  a  variable  medium 

(13) 

according  to  Fermat's  principle  is,'  ' 


d  / 


^  •  n  =0 


(2.  34) 


where  Is  a  unit  vector  along  the  ray,  Is  the  distance  along  the  ray  and  n  is 


the  index  of  refraction.  Equation  (2.34)  can  be  rearranged  to  read, 


log  (7^)J  -T{T‘  5  log  . 


(2.  35) 


To  evaluate  these  operations,  Introduce  the  decompositions: 


(2.36) 


where  "t  Is  the  initial  value  of  "ttyi  is  the  mean  index  of  refraction  and  is  very 
o  .  (  o 

TT_ - 

'  'S.  Chandrasekhar ,  Monthly  Notices,  Royal  Astro.  Soc.,  112,  475  (1952). 

See  also  H,  W.  Liepmann,  "Reflection  and  Diffusion  of  a  Light  Ray  Passing 
Through  a  Boundary  Layer,  "  Douglas  Aircraft  Co.,  Inc.,  Rpt  SM- 14397,  May 
16,  1952,  for  a  discussion  of  a  similar  phenomenon, 

Kerr,  "Propagation  of  Short  Radio  Waves,  "  Rad.  Lab.  Series,  Vol.  13, 
McGraw-Hill,  1951,  page  44. 
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-carlv  equHl  to  unity,  and  N,  3.8  the  variable  part  ofri.*  To  first  order  in  e  and 
N'j  (which  must  be  very  smaU),  Eq.  (2.35)  can  be  written, 


^  '  ' o - 

‘  °  <:)J: 


(2.  37) 


de  I  djl  defines  a  differential  angle  (and  direction)  of  ray  displacement.  Equation 
(2.  37)  g  ves  i.his  angle  in  terms  of  the  gradient  normal  to  the  path;  integrating 
over  the  whole  pa  t^-  gives  the  total  angular  deflection. 

Let  us  return  to  the  rotating  line -of-sight  problem.  Using  cylindrical 
coordinates  (with  the  z-axis  normal  to  the  page),  the  orthogonal  components  of 


angle  dfviation  obr-iined  from  Eq.  (2.  37)  are: 


^  f  dj?  - 1- 

*"0  d  z 


E  1  c5N, 


Vq  =  d;:-- — L  . 

>-''0  SQ 

The  auto-correlations  of  these  two  functions  are  thus: 


+  t)>  =  d/j  f  dj?2 


\  3z  ,  3  z 


(2.  38) 


(2.  39) 


(2.40) 


<;*  (t)«  (t  +  T)>  r^dA— i— 


(2.41) 


Considering  Eq.(2.40t  first,  we  have  by  differentiating  inside  the  average  signs 


in  the  definition  of  p  . 


/-  d  N ,  I  d  N ,  \  , 

/  1  I  I  X  _  /.,2  \ 


^  -  <N^,: 


rj^  Pi^iz) 


(2.42) 


^  In  terms  of  N  unat.s,  N  j  -  AN-  10 
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with  rj2  given  by  Eq.  (2. 43)  below  and  ^  equal  to  •  In  cylindrical 

coordinates,  the  distance  between  two  general  points  (1)  and  (2)  is: 


12 


1/2 


(2.43) 


Evaluating  Eq.(2.42)  in  the  x,  y  plane  (zj  =  =  0),  one  obtains 


'aN, 


9  z 


aN, 


a  z ! 


z.=0 

z^=0 


=  -  <Ni> 

"12 


dp(rj2) 


(2.44) 


12 


Using  Eq.(2.2)  and  taking  0j  =  0,  we  have 


<^^(t)^z(^  +  t)>  = 


^1  ‘  ^^1^2 


X  S 


1/2 


■^o 


^^Z  *  ^V2 


1/2 


(2.45) 


The  methods  used  in  integrating  Eq.  (2. 21)  may  be  applied  here  to  yield 


(t)^  (t  +  t)>=  tt<2^n'^>  .  10^^  . 


z'  ' 


4 


+  K_j(^)L^(^)  (2  46) 
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With 


^  = 


L9|t 

eo 


It  will  be  noted  that  Eq,  (2.46)  becomes  infinite  as  t  — ►  0;  which  is 

to  say,  the  mean  square  deviation  of  is  infinite.  This  is  simply  a  conse¬ 
quence  of  our  choice  for  the  correlation  function  made  in  Eq.  (2  2).  Unfortu¬ 
nately,  Eq.(2.2)  describes  a  (Markov)  process  which  has  no  derivatives.  This 
need  not  disturb  one  unduly,  since  a  spectrum  can  still  be  found  from  Eq.  (2  46) 
by  taking  its  cosine  transform.  The  spectrum  W^(f)  obtained  from  the  cosine 
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Figure  8 

Angle  of  Arrival  Spectra  for  Rotating  Line -of -Sight 


transform  of  Eq.  (2. 46)  is 


Ztt  •  10" 


W^(f)  = 


e 


log 


W  , 


(2.47) 


with  Wj  defined  as  in  Eq.  (2.  31).  This  is  plotted  in  Fig.  8.  This  spectrum  is 
meaningful,  except  as  the  frequency  f  oo.  The  behavior  at  large  f  is  such 

pOO 

that  df  W{f)  diverges,  yet  the  density  at  low  frequencies  is  correct.  Since 

the  data  arc  too  crude  to  indicate  a  better  choice  for  p(r),  we  have  tried  to 
extract  the  maximum  amount  of  usable  information  from  Eq.(2.2).  A  different 
choice  for  p(r)  might  well  give  a  spectrum  whose  integral  converges,  yet  the 
low  frequency  density  would  not  be  more  trustworthy  than  that  derived  from 
Eq.  (2.47).  One  might  thus  be  led  to  a  false  sense  of  security  in  calculating 
the  total  mean  square  as  a  definite  number;  when,  in  fact,  slight  variations  in 
p(r),  which  might  fit  the  data  equally  well,  could  produce  large  changes  in 


The  auto-correlation  of  Pg  and  the  corresponding  spectrum  can  be 
found  by  a  procedure  similar  to  that  above. 


<Pg(t)pQ(t  +  T)  > 
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^<AN^> 
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(2.48) 


and 


Wg(f)  = 


2tt^  •  10'^^<AN^"^  1 


(0)' 


2TTf  A 
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0L  / 


'  J 
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Equation  (2. 49)  is  also  plotted  in  Fig.  8, 


(2.49) 
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C  Probnbililv  Dcrsit^e s 


(1)  Total  Fiold  Strength 

( I  ‘ 

In  ‘he  der  vat:on  '  of  Eq.  (2. 3>  it  is  assumed  thjjt  the  volume  elements 

to  which  the  equation  applies  are  ’arge  in  linear  dimensions  compared  to 

When  this  is  true,  the  scattering  from  separate  volume  elements  will  have 

random  phasing  Tie  total  field  is  the  sum  of  t^'ese  scattered  fields  and  its 

amplitude  wjl  be  Rayleigh  distributed.  The  direct  ray  has  random  phase  but 

constant  ampl'^ude.  The  sum  of  direct  and  scattered  waves  is  then  the  sum  of 

(14)(15) 

a  sine  wave  plus  Gaussian  noise  and  has  the  well  known  form, 

E  i  /  EP  \ 

p(E)  =  —  expl - I  f  -  ).  (2.50) 

Q  .  20  J  °  V  Q  / 


Here  p{E)  is  the  probability  density  of  finding  the  sum  amplitude  E  between  E 
and  E  +  dE,  Q  is  the  total  scattered  power,  P  is  the  direct  wave  amplitude, 
and  Iq(x)  is  a  modified  Bessel  function  of  the  first  kind  The  case  P^^  Q  . 
which  is  of  interest  here,  allows  Eq.(2.50)  to  be  reduced  to 


p'E)  ^ expF.iE  -  P)^/2Q  . 

v/TttQ  L 


(2.51) 


This  shows  that  the  received  amplitude  will  be  very  nearly  normally  distributed, 
with  mean  P  and  standard  deviation 

(2)  Phase  Angle  Probability  Densities 

As  just  stated  above,  the  phases  of  the  scattered  waves  will  be  random, 
i.e.,  Q  is  uniformly  distributed  from  -tt  to  tt.  Th  s  must  be  combined  with  the 
phase  deviation  of  the  direct  ray.  From  Eq.(2.  17),  the  total  phase  delay  is 


^Rice,  Mathematical  Ana.lysis  of  Random  Noise,  BSTJ,  £4,  46  (1  945). 
(15) 

'  'Lawson  and  Uhleiibeck,  Threshold  Signals,  Rad.  Lab.  Series,  24, 
McGraw-Hill  (1950). 
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given  as  the  sum  of  the  individual  phase  delays  in  the  turbulence  clusters.  The 
sum  will  be  normally  distributed  by  virtue  of  the  central  limit  theorem. 

This  theorem  predicts  a  normal  density  function  over  an  infinite  range  for  phase, 
in  the  more  usual  case  where  phase  is  considered  modulo  Zrr,  the  resulting 
density  function  is'  ' 

p(a)  =  1.  0  (2.  52) 

Zn  Vti  Ztt  / 

where  0^  is  a  theta -function.  When  <<C^  Ztt  ,  this  reduces  to  the 

normal  form 

2 

g 

p(q)  =  -_! -  e  ^  .  (2.53) 

v/ztt  'Cq^> 

The  angle  of  arrival  density  for  the  scattered  wave  is  given  by  Eq. 

(2.  13).  The  ray  deviation  density  will  be  normal  in  two  dimensions  (and 
uncorrelated  in  those  two  dimensions). 

D.  Numerical  Estimates 

Some  estimates  have  been  made  for  the  scattering  parameters  which 

(18) 

appear  in  the  foregoing  expressions.  Very  near  the  ground,'  '  one  finds 
AN  i  10.  Microwave  refra ctometer  flights^^^  indicate  that  AN  varies  between 
1/2  and  IN  units  from  several  hundred  to  thirty-thousand  feet.  This  altitude 
range  includes  most  of  the  tropospheric  region  in  which  refractive  fluctuations 
affect  propagation.  The  correlation  length  has  a  ground  level  value  of  20 

Cramer,  "Mathematical  Methods  of  Statistics,  ”  Princeton  University 
Press,  1946. 

(171 

'  'P.  Levy,  Societe  Math,  de  France,  Comptes  Rendus,  p.  32,  1938. 

M.  Crain  and  J.R.Gerhardt,  Proc.  I.R.E.,  40,  50(1952). 
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feet  and  increases  rapidly  to  an  asymptotic  value  of  one  or  two  hundred  feet. 

Since  the  preciSf  values  of  and  AN  vary  with  both  geographic  location  and 

time,  we  shall  choose  the  convenient  averages  AN  =  1/2  and  X  =100  fe*  t. 

o 

We  take  our  line -of-sight  transmission  path  length  L  to  be  10  000  feet  and 
assume  an  antenna  beamwidth  (3  of  3°  We  shall  also  take  \  =  6  cm  so  as  to 
(later)  validate  the  Raylcigh  approximation  and  avoid  tie  attenuation  problems 
characteristic  of  shorter  wavelengths. 

Contributions  to  the  RMS  phase  error  arise  from  both  the  scattered 
and  direct  wave.  Substituting  the  above  parameters  into  Eq.  (2.  1  1)  and  Eq 
(2.19).  we  find- 

2 


and 


Q  ?  =  2.5  •  10  ^  radians^  . 

sc 


Q,  =  0.  5  •  10  ^  radians^  . 

dir  ect/ 


(2  54) 


Scintillation  of  the  direct  wave  is  evidently  the  dominant  effect  here  This 
error  represerMs  one-tenth  of  a  cycle  at  5,000  Mes  thereby  posing  a  genuine 
problem  for  systems  requiring  high  p.base  stability  over  the  transmission  path. 
The  spectrum  of  such  phase  fluctuations  is  of  considerable  interest 
also.  S>nce  the  characteristic  turbulent  self-motions  (and  hence  doppler 
shifts)  are  known  only  in  terms  of  the  qualitative  picture  of  eddy  degeneration 
(Refs.  (3)  and  (4))  we  have  computed  spectra  for  a  moving  line -of-sight .  We 
have  just  considered  a  rotating  path,  the  spectrum  of  which  is  given  in  Figure 
7.  A  measure  of  spectral  width  is  given  by  2iTf/^^/L0  =1.  If  one  assumes  that 
the  line  -  of- sight  is  rotating  so  as  to  track  an  aircraft  flying  1,000  fps  at  an 
altitude  of  10,000  feet  (overhead),  0  is  approximately  0  1  radians/soc.  Using 
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this  with  the  parameters  above,  the  effective  width  is  1.59  cps. 

A  somewhat  different  model  is  obtained  by  anticipating  the  result  of 
Part  IV  which  gives  the  spectrum  for  a  line -of-sight  moving  parallel  to  itself 
through  a  stationary  turbulence  configuration  with  speed  v.  ♦ 


One  may  identify  v  with  a  local  wind  speed,  which  blows  the  (stable)  turbulence 
structure  through  the  line -of-sight.  This  expression  is  halved  when  f  =  0.  122 
v//{f^,  so  that  with  v  =  20  ft/ sec  and  =  100  feet,  an  effective  (noise)  band¬ 
width  B  =  0.  02  cps  is  produced.  Very  near  the  ground,  one  has  \0  feet, 

so  that  B  0.2  cps. 

The  received  power  is  also  smeared  over  a  small  solid  angle  centered 

on  the  line-of-sight.  Off-axis  turbulent  clusters  scatter  energy  into  the 

antenna  according  to  the  angular  distribution,  Eq.  (2.  13).  This  ratio  is  reduced 

\  1-3 

by  a  factor  of  four  when  0  becomes  -  —  10  radians,  which  represents  a 

21.4  3 

very  narrow  pencil  of  rays  indeed. 

We  have  also  investigated  the  angle -of-arrival  fluctuations  which  the 
direct  ray  experiences.  The  exponential  correlation,  Eq.(2.2)  produces  an 
infinite  RMS  deviation  because  of  its  cusp-like  behavior  at  the  origin.  In 
general,  the  spectrum  widths  for  deviation  of  the  rays  is  of  the  same  order  of 
magnitude  as  that  for  the  phase  variations.  It  would  be  of  great  interest  to 
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repeat  the  foregoiniz  calculat.ons  when  a  better  correlation  function  is  established 
on  physical  grounds. 


III.  RAINFALL 

The  theory  of  mui*  path  transmission  through  singly  scattering  media 
developed  in  Pait  IJ  maybe  exploited  to  discuss  rain-produced  scintillation 
errors  as  a  radar  beam  traverses  a  precipitation  area  The  drops  w  11  be 
idealized  as  perfect  dielectric  spheres.  The  entire  signal  perturbation  thus 
results  from  c  ass. raj  scattervig  and  no  phase  d  sturb.ince  of  the  direct  wave 
is  considered.  When  these  results  are  combined  with  appropriate  m<  ’corologi 
cal  data,  the  theory  predicts  bo*h  angular  distributions  and  phase  error 
frequency  spectra  for  the  received  power. 

A.  Static  Mult -path  Effects 

Ou’  descr-pt'on  begins  with  an  expression  for  the  scattering  by  a 
single  drople*,  as  lllastr  itcd  in  Fig.  9.  Since  the  raindrops  are  considered 
spherical  one  m\»5t  deduce  tic  r  itio  of  scattered-to-  ncidcrt  power  for  a 
sphere  of  radius  a  and  dioioetr.c  constant  e  .  Th.s  qaaiit.lv  has  been  computed 
by  classical  electromagnetic  theory  in  many  places  Wp  shall  see 

presently  that  most  ra.ndrops  are  less  than  2  millimeters  in  diameter,  and, 
since  we  rons.der  microwave  propagation,  one  may  pass  directly  to  the 
Rayleigh  limit  ^Ztra  ' \ <C3C  1 '  ^ 

7T9T 

'  'D.E.Kerr,  Propagation  of  Short  Radio  Waves,  MIT  Rad.  Lab.  Series, 

New  York,  McGraw-Hiil,  Jl_3,  445  (1951). 

.  A.  Stratton  Electromagnetic  Theorv,  New  York,  McGraw-Hill,  p.  563 
(1941). 


?0  - 


(3.  1) 


S(R)  _ 
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4tt  •  R 
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^  /  2TTa 

•  ^ 
sin^(0)  +  cos^(6)  cos^(0) 

je  +  2 

representing  reradiation  by  the  drop  as  an  electric  dipole.  The  scattering 
cross  section  per  droplet  per  unit  solid  angle  is  obtained  by  multiplying  this 
ratio  by  4TrR^.  To  compute  the  scattering  cross  section  per  unit  volume, 
which  appears  in  the  multipath  density  £q.(2.  8),  one  must  estimate  the  number 
of  drops  per  unit  volume  in  each  radius  interval  (da).  Rigby  and  Marshall  have 
fitted  the  exponential  function  £q.(3.2)  to  meteorological  experiments .  ^ ^ 


dN 

dadv 


c 


drops 
(cm)(cm  ) 


(3.2) 


Figure  9 

Scattering  of  a  Plane  Wave  by  a  Spherical  Drop 
Their  normalization  constant  c  is  roughly  0.08  and  the  "drop  scale  size 
parameter,  "  b,  is  given  in  terms  of  the  rainfall  rate  Pv  (measured  in  milli¬ 
meters  per  hour),  by 


^^£.  C.  Rigby  and  J.  S.  Marshall,  "Modification  of  Rain  with  Distance  Fallen,  " 
Scientific  Rpt.  MW-3,  MacDonald  Phys.  Lab.,  McGill  University,  Jan.  1952. 
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R 


0.21 


(cm) 
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(3.3) 


Table  1  correlates  R  and  b  with  familiar  precipitation  intensities. 


Table  I 

Rainfall  Rate  Versus  Drop  Scale  Size 


Rain  Type 

R 

b 

Drizzle 

0.25  mm/ hour 

0.018  cm 

Light  Rain 

1.00 

0.024  " 

Moderate  Ra^n 

4.00 

0.033  " 

Heavy  Rain 

10.00 

0.040  " 

Very  Heavy  Rain 

16.00 

0.044  " 

It 

20.00 

0.046  " 

tt 

25.00  " 

0.048  " 

1 1 

30.00  " 

0.050  " 

It 

50.00 

0.055  " 

It 

100.00 

0.064  " 

tt 

150.00 

0.070  " 

Expression  (2.8)  for  the  static  multipath  power  distribution  may  now 
be  computed  by  combining  Eq.(3.  1)  and  (3.2)  and  averaging  over  all  drop  sizes. 

J'" da  c  jde  sinS  ('^"d?  4-ra^  ^ 

P  o  o 


e  +2  \  \  J 


1^8  in 


^(0)  +  COS^(0)  COS^(0) 


X  -  r  (1  -  cos6) 


]■ 


(3.4) 


A  spherical  coordinate  system  erected  on  the  receiver's  pointing  direction  has 
been  used  to  write  out  the  beam-limited  (primed)  volume  integral  (see  Figure 
10).  For  a  narrow  beam,  one  may  disregard  the  entrance  disparities  between 
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Figure  10 

Coordinates  for  Scattering  by  Rainfall 


( 1 )  Phase  Delay 

With  Eq.  (3.4)  and  the  superposition  of  Figure  3,  one  computes  the  RMS 
phase  error  in  a  sinusoidal  voltage  received  at  the  origin  as: 


:a^>  =  C  dx  sin^  f  .  QM  . 

o  \  \  /  P 

o 

r*^d0  sinO  [l  +  co8^{0) 

\  \  J  L  J 


c  -  1 


e  +  2 


(continued  on  following  page) 
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■  c  o  s  (  0  ^ 


(Sec  0-1)  sin 


(Sece  -  1) 


4Tr/  X. 


4Tr/X. 


(3.  5) 


The  trigonometric  terms  may  be  ignored  in  comparison  wiCh  the  viewing  thick¬ 
ness  (R^  -  Rj);  and  the  dielectric  constant  (e)  assumed  constant  to  our 
approximation. 


<a^>  .  (R^  -  Rj)  ^  (P)  . 


(3.6) 


^  denotes  a  beam  shape  factor, 


-f /Q\  'i®  sin(e)  r,  2,q. 

3  (p)  =  (  - —  1  +  cos  (9)  , 

^o  cos(0)  L 


^P^  for  P«1  .  (3.7) 

Our  result  indicates  that  varies  about  as  the  seven-tenths  power  of  the 

rainfall  rate  and  directly  with  the  square  root  of  the  total  number  of  participating 
droplets.  This  latter  dependence  is  characteristic  of  incoherent  scattering 
processes  and  might  have  been  anticipated  from  the  original  assumptions.  To 
provide  numerical  estimates  for  Eq.(3.6),  we  note  that  the  dielectric  constant 
of  water  varies  from  e  =  78.5  -  il2.3atX  =  10  cm  to  34.2  -  i35.9atX  =  1.24 
cm  for  t  -  18°C.^^^^  The  corre  spending  values  of  I— -t-tI  range  from  0.  9286 
to  -A  .  9206,  and  we  shall  choose  the  average  value  0.  925.  In  Figure  1 1 ,  vve  plot 
■  (R^  -  R^)  for  various  wavelengths  (X)  and  rainfall  rates  (R).  From 
those  curves  we  read  for  X  =  6  cm  and  R  =  100  mm/hour  (i.e.,  a  very  heavy 
rain) , 


Kerr,  op.  cit.  ,  page  610. 
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(3.8) 


<.a^  "  ' :3.4.  5  •  1  0'^  (m‘ P^(R2  -  R ^ )  . 


A  rvpica]  antenna  beam  p  =  4°  and  rainstorm,  R^  -  Rj  =  1000  meters  predicts, 


°RMS  —  ^‘^3  •  10  ^  radians  . 


(3.8a) 

This  error  is  an  order  of  magnitude  larger  than  the  scattering  effect  of  static 
tropospheric  turbulence  derived  in  Part  II.  The  sensitivity  of  Eq.  (3. 6)  to 

rainfall  rate  (R)  encourage  one  to  estimate  local  meteorologi¬ 
cal  'conditions  with  some  care  before  attempting  precision  measurements. 

( ^  '  -*^-nguIar  Distribution 

An  angular  distribution  of  power  about  the  receiving  antenna's  pointing 
dirrction  is  predicted  by  the  foregoing  theory.  If  one  integrates  expression 
(1.  '  i  over  a)i  path  differences  (x) ,  the  total  scattered  power  becomes: 


p 

d0  sin(0 

oo  -)  /u  R,Sec0 

)  r  r  ^  dr 

2rr  j 

S 

e  -  1 

o 

o  R^SecO 

o 

e  +  2 

Ztra 

K 


2  2  2 
sin  (0)  +  COB  (0)  cos  (0) 


(3.9) 


The  angular  dependence  of  this  scattered  power  is  recognized  as  the  integrand 

of  the  cos{0)  integration. 


dP 


417^(6);  /Zirb 
- - -  c 


P^  d(cos0) 


b  (R^-R,) 


1  +  coa  ley 
L  cob(0) 


(3.  10) 


The  angle -dependent  factor  is  plotted  in  Figure  12,  where  a  relatively  flat 
distribution  for  small  beamwidths  is  exhibited.  Increasing  values  of  this 
function  for  large  6  are  due  to  the  inclusion  of  an  ever-increasing  number  of 
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scatterera — the  number  going  to  infinity  as  near  horizontal  viewing  is  included. 
A  realistic  consideration  of  attenuation  effects  rectifies  this  problenri.  We  shall 
continue  to  use  the  natural  cut-off  of  narrow  beam  antennae. 


Figure  12 

Angular  Distribution  of  Rainfall-Scattered  Power 

B.  Scintillation  Spectra 

In  the  foregoing  we  derived  an  expression  for  the  RMS  value  of  the 
multipath  phase  error.  Of  equal  importance  is  the  frequency  spectrum  of  the 
phase  scintillations.  The  anomolous  signal  is  compounded  from  the  scatterings 
by  thousands  of  raindrops  in  the  beam.  Time  variations  of  the  phase  error 
may  be  regarded  as  the  composition  of  Doppler  frequency  shifts  from  the 
individual  droplets  moving  relative  to  one  another  in  changing  order.  We  may 
anticipate  the  magnitude  of  these  frequencies  from  the  Doppler  relation: 
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Li  -  u!  c  .  For  relative  drop  speeds  of  one  meter  per  second  and  f^  =  5000 
mcs,  Li  -  16  cps  and  such  a  scintillation  is  evidently  of  some  importance  in 


precision  measurements. 

We  first  resolve  all  drop  velocities  along  the  pointing  (z)  direction  of 
the  antenna,  as  shown  in  Figure  12.  If  the  incident  plane  wave  meets  a  drop 
at  dv  with  speed  u,  the  drop's  motion  away  from  the  incident  wave  lowers  the 
incident  frequency  (seen  by  the  moving  drop)  by  an  amount  u/ c  f^  .  The  drop 
reradiates  thio  shifted  frequency  at  an  angle  0  to  the  incident  beam  thereby 
causing  an  increase  f^  u/ c  cos(6).  The  net  frequency  change  during  the 
scattering  process  is^f  =  f^  •  u/c  •  (1  -  cos0).  We  shall  construct  a  joint 
path  difference-frequency  power  spectrum  for  the  multipath  signal  scattered  by 
the  volume  element  dv  at  (r,0,0)  moving  with  speed  u  parallel  to  the  incident 
beam.  By  using  the  analytical  device  introduced  in  Eq.  (2.6)  of  Part  II,  one 
may  show  that  the  ratio  of  the  scattered -to-incident  power  with  path  difference 

between  x  and  x  +  dx  and  signal  frequency  between  f  and  f  +  df  is  given  by. 

00  00  -I  ,  r  "I 

-  '  =  ^  da  N{a)  ^  du  P{u|a)  J  ~  <y;(a,r,0)  6 |x  -  r ( 1  -  cos0)j 


f  -  f  -  (1  -  cos 
o  o  c 


e)j  . 


(3.  11) 
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Figure  13 

Drop  Velocity  Resolution  for  Doppler  Shift  Analysis 


A  prinrte  on  the  volume  integration  restricts  that  process  to  the  radar  beam. 
P(uja)  denotes  the  conditional  probability  distribution  for  a  velocity  component 
u  parallel  to  the  auteniia  axis  for  drops  with  radius  a.  We  do  not  know  this 
density  a  priori,  and  its  determination  is  the  central  problem  of  this  section. 

We  shall  see  presently  how  one  may  infer  P(ula)  from  other  propagation 
experiments — notably,  backscattered  (echo)  fluctuations  from  rainfall. 

The  integrations  of  Eq.  (2.  54)  may  be  performed  with  relative  ease 
if  one  assumes  P(u'a)  to  be  independent  of  drop  size  (a).  Such  an  assumption 
is  probably  warranted  if  the  relative  drop  motion  is  due  to  local  gusts  and 
turbulent  conditions  in  the  rainfall  area.*  In  a  C.  W.  phase  comparison  system, 
the  scintillation  frequencies  are  given  by  excursions  from  the  carrier  frequency 

It  is  also  possible  to  obtain  a  Doppler  shift  from  the  differential  speeds  with 
which  drops  of  different  sizes  fall.  This  effect  can  be  calculated  from  Stoke 's 
law  for  the  fall  of  spheres  in  a  viscous  medium.  We  shall  return  to  thir  caie 
later  for  comparison. 
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,  80  that  the  frequency  spectrum  is  a  function  of  v  =  f  -  f^  .  Inserting  the 
appropriate  forms  for  N(a)  and  <T  from  Eq.  (3.  2)  and  Eq.  (3.  1)  respectively,  one 
finds . 


^  /  1  \  P  ^  *>  2  TT 

— ^  du  P(u)  ^  da  c  e'  ^  ^  d6  8in{6)^J^  dr  ^  d^ 


P  o 

o 


Rj  o 


r  2  4  -| 

•<|^4Tra^  ^  ^  I  ^  {8in^(0)  +  co8^(^)  C08^(e))  > 


•  6  X  -  r(l  -  cosO)  •  6  ^  ~  0  "  coa{Q) 


(1)  Phase  Dela' 


(3.  12) 


If  one  uses  this  expression  to  compute  the  phase  error,  the  frequency 
spectra  for  emerges  without  further  effort. 


<a^(v)>  =  J^°°dx  sin^ 

o  '  X.  / 


Q(x,y) 


X  /  P 
o 


-R,)g(>^)~  rde  -  1  +  C08^(e) 

2  ^  f  1  -  co8(e)  L 


•  P  - - 1  .  (3.13) 

_f^(i  -  co8(e) 

The  frequency  dependence  of  the  forward  scattering  is  apparent  from  Eq.(3.  13). 
Since  the  velocity  distribution's  argument  c/f(l  -  co80)  ^  is  large  for  the  small 
values  of  (1  -  cosQ)  inherent  in  a  narrow  beam,  the  net  effect  is  to  relegate 
most  of  the  scintillation  power  to  small  frequencies.  g(X  )  denotes  a  wave 
length-dependent  term  arising  from  the  drop  radius  averaged  cross  section. 


40  - 


g(X)  = 


da(ct"^^^^)  4na^ 


-  1  /  Zffa  \^ 

€  +  21  \  \.  / 


=  0.08  4irb^  LIJlJi  • 

\  X.  /  le  +  2I 


(6)! 

27 


ib^ 


(ft 


(3.14) 


We  assume 


ie  -  1 


!€  +  2l 


^  0.925,  as  before,  and  b  is  given  in  terms  of  the  rainfall 


rate  R  by  Eq.  (3.  3). 

To  determine  the  function  P(z),  we  shall  first  exploit  experimental 
data  on  spectral  analysis  of  precipitation  echoes. 


Figure  14 

Echo  Process  for  Moving  Drops 


Figure  14  exhibits  the  ba ckscaiter :ng  process  for  a  moving  drop,  A  plane  wave 

sent  out  by  the  antenna  0  radians  off-axis  enjoys  a  frequency  shift  u/ c  cosG  , 

both  on  reception  and  reradiat  on.  The  appropriate  variation  of  Eq.  (3.  1  1)  for 

the  backscat^er cd  ampl.tudes  is 
00  00 


Jda  N(a)  J’duPCulajj’ 


P 

o 


X  -  r  (i  -  cosG) 


J  •  6  -  2f^  ^  cosG  J 


(3.  15) 


An  amplitude  scint  llation  spectrum  may  be  deduced  from  this  form  by  inte¬ 
grating  over  (x)  and  substituting  the  appropriate  backscatter ing  cross  section 
from  Eq.  (3.  1). 


^  y  da  N{a)4TTa^  | 
^  °  o  ! 


e  -  1 


€  +  2 


2TTa 

X 

oo 


4  ^oo 


■)  ^ 


du  P{u) 


P  R,  2Tr 

C  d0  sjn(O)^^  dr  ^  d^  dx  6 
R  j  o  o 


X  -  r{l  -  cos 


6^ 


6  X/  -  2f 


o  7 


Pog(>‘)’'(R2  -  Rj)  ^  J’ de  un(e)  p 


C  7' 


2f  cos(9) 
-  o  '  ' 


The  spectrum  normalized  to  unity  at  i-'  =  0  is  simply, 

P 

dG  ta n(0)  P 


G{t»  = 


C  7' 


2f  cos(G) 
u  o 


P(o)  ^n  ,cos{p) 


'J 


(3. 16) 


(3.  16') 


The  small  ( <4  )  values  assumed  for  p  permit  one  to  replace  cosG  by  unity  in 
Eq.(2.6l)  and  thus  extract  P(z)  from  within  the  integrand. 
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(3.17) 


This  important  relation  permits  one  to  infer  the  velocity  distribution  of  rain¬ 
drops  from  the  experimental*  e.ho  fluctuation  spectra  G(t.).  In  Figure  15  we 

(19) 

reproduce  the  measurements  of  Goldstein  and  others  for  X  =  9.2  cm,  as 
measured  on  three  occasions.  These  cu:  ves  give  just  the  normalized  (Doppler) 
frequency  power  density  of  the  backscattered  amplitudes,  G(^). 


Figure  15 

Power  Frequency  Spectrum  of  the  Fluctuations  of 
Precipitation  Echo  on  9.2  cm  as  Measured  on  Three  Occasions 


Note  that  Eq.(3.  17)  is  independent  of  the  beamwidth  (P)  used  in  these 
experiments . 

E.  Kerr,  op.  cit.,  page  576. 
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Our  integral  expresaion  (3.  13)  for  the  phaae  error  may  now  be  combined  with 
Eqa.  (3.  16')  and  (3.  17)  to  yield: 

'  =  I.(R,  -  R  )g(v.)  . -S- P(o)  j-P  lL£}il(£L. 

'  2  .  f  ‘-'o  I  -  coi(9) 

o 

•  fl  +  coa^{0)l  •  c[—^ - 1  .  (3.  18) 

I-  J 

The  remaining  integration  haa  been  performed  numerically  uaing  the  valuea  for 
G{  i/)  giv'cn  by  Figure  15.  In  Figure  16,  we  plot  the  normalized*  apectrum  for 
aeveral  beamwidtha  (p).  Theae  curvea  emphaaize  the  concentration  of  power  at 
low  frequenciea. 

One  may  alao  compute  a  phaae  angle  apectrum  by  ( ona.aering  the 
detailed  motion  of  individual  dropa.  Imagine  an  enaemble  of  dropa  falling 
through  atagnant  air  under  gravity.  The  viacoua  drag  force  predicta  greater 
fall  apeeda  for  heavy  drops  than  for  light  onea.  Extremely  amall  spheres 
(a  <  0.  008  cm)  obey  Stoke 'a  law  for  ateady  state  terminal  velocity,  ♦♦ 


8 

=  _  g  a 
3 


-  1 


'D 


Very  large  drops  are  relatively  unstable  and  tend  to  break  up  into  smaller  ones 
during  long  descents.  The  intermediate  region  is  of  some  interest  for  this 
study,  in  view  of  the  drop  scale  sizes  listed  in  Table  I.  Fortunately  enough. 


Normalized  so  that  d  »  \q  j^)  =  1 . 

o 

♦  ♦ 

Cj^  denotes  the  drag  coefficient,  the  air  density,  the  drop  density;  a  is 
the  drop  radius  and  g  the  acceleration  of  gravity. 
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'22; 

Gunn  and  Kinzer’  have  measured  the  terminal  velocities  of  drops  with  diame¬ 
ters  in  the  range  0.  01  <  D  ^0.  68  cm.  Their  results  may  be  summarized  by  an 
approximate  formula, 


v(a} 


-  /ia 
e 


(3  19) 


with  K  =  1020  m/sec  and  A  -  5.04  cm”^.  The  conditional  probability  density 
P(u|a)  for  the  distribution  of  velocities  with  fixed  radius  (a)  may  be  given 
explicit  form  in  terms  of  Eq.  (2.  19). 


P^uia)  =  ^  *  v(a)J 

If  the  angle  between  the  local  vertical  (i.e.,  it  )  and  the  antenna's  pointing 
direction  is  denoted  by  7,  then  u  =  v  cos(y).  Equations  (3.  13),  (3.  19),  and 
(3.20)  may  then  be  combined  to  give  the  phase  angle  spectrum: 


(3.20) 


<'aMy  =  J’  d 


.  Z/ZlTX 

X  sin  (  - 
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.00  00  r 

daN(a)  du  6  u  scc(y)  -  v(a) 
/  o  o  L 
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Rz 


•  p  do  8-n(0)^J^  ^  cr(0,0;a)  6 


r 


R  j  o 


X  -  r(l  -  COS0) 


) 
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7'  -  f  —  (1  -  COS0) 

O  C  '  'j 


1 


2:4tt^(0.08)(R,  -  R.) 

^  ^  f  ie  +  2 

o 


2  00 


S 


^-(2a/b) 

v(a) 


^^  R.  Gunn  and  G.  D  Kinzer,  Jour,  of  Meteorology,  243  (1949).  Their  ex¬ 
periments  were  performed  under  760  mm  pressure,  with  a  temperature  of 
20®C  and  relative  humidity  of  50  per  cent. 
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The  scintillation  frequency  (iJ)  persists  in  the  (drop-size)  integration  condition, 


v(a)>  SJL  ■«<=(Y>  .  (3.23) 

f^  1  -  C08(p) 

If  the  curve -fitted  function  (3.  19)  is  inserted,  the  (origin-normalized)  phase 
spectrum  becomes: 


F(v)  = 


J’ 


00  6  -(2a/b) 

da  ^ - 
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~  A  a. 


00  6  -(2a/b) 

a  e 


J' 


da 


1  -  e 


-  /A  a 


(3.24) 


with 


^  =il„  1  .  _£££(ir]. 

L  Kf^  (1  -  cospy 


It  is  clear  from  this  form  that  all  frequencies  greater  than 


- —  f  cos(v)(l  -  cosB) 
max  o 

c 

arc  eliminated  in  Eq.  (3.  24).  We  may  estimate  this  cut-off  by  assuming: 
f  =  5000  mcs,  cosv  =  0.  4,  p  =  4°,  K  =  10  m/sec,  so  that  v  =16  cps. 

The  quotient  (3.24)  was  computed  numerically  for  moderately  heavy  rain, 

R  =  4  mm/hour  and  is  displayed  in  Figure  17  for  the  range  0.  01<  <  16  cps. 

Rather  less  realism  may  be  associated  with  this  model  than  the 
radar-precipitation  echo  approach.  This  statement  is  based  on  the  common 
experience  that  heavy  rainstorms  are  fraught  with  turbulent  gusts,  each 
probably  sufficient  to  dominate  the  laboratory  result  (4.  19).  It  would  be  very 
valuable,  however,  to  perform  a  controlled  experiment  measuring  line -of -sight 
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&UX-  ocAiiering  Dy  ir'alling  Drops 


PC  Sec  y 


and  echo  scintillations  simultaneously  within  the  same  rainstorm.  Such  a  test 
could  criticize  the  result  (3.  16')  and  perhaps  give  information  on  the  mixture  of 
turbulent  and  free-fall  drop  motion  within  rainstorms. 

C.  Ray  Bending  in  Rainstorms 

In  addition  to  scattering  off-axis  energy  into  the  antenna,  a  rainstorm 
will  act  as  a  prism  and  bend  rays  passing  completely  through  it.  If  the  drop 
sizes  are  small  compared  with  the  wavelength  of  the  incident  radiation,  the 
equivalent  index  of  refraction  for  a  collection  of  perfect  dielectric  spheres 

I €  1  i  ®  • 

=  1  +  Zrra^N  .  (3.25) 

Here  a  is  the  sphere  radius  and  N  the  number  of  spheres  per  unit  volume.  One 
may  use  Eq.(3.2)  to  average  this  expression  over  all  drop  sizes. 

^  0.  2b^  (3.26) 

For  a  heavy  rain  R  =  10  mm/hour  ,  b  =  0.04  and 

^  ^  1  +  1  .  10'^  .  (3.27) 

The  bending  at  entrance  and  exit  of  a  ray  passing  through  such  a  region  can  be 
computed  with  standard  techniques. 

When  the  drop  spacing  is  small  compared  with  the  wavelength,  one  must 

(23) 

exploit  a  detailed  picture  of  the  scattering  processes.  H.C.  Corben  has 
performed  such  an  analysis  and  concludes  that  the  essential  factors  in  Eq.(3.25) 
are  preserved.  He  finds  that  multiple  scattering  may  be  disregarded  for  most 
applications . 

TZD - r~ 

'  'Private  communication,  August,  1954. 
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IV.  CLOUDS 


Index  of  refraction  data  for  clouds  are  notably  sparse.  From  meteor- 

logical  measurements  and  aircraft  gust  load  experience,  However,  it  is  known 

that  the  interiors  of  cloud  formations  are  very  turbulent  and  contain  relatively 

large  index  fluctuations.  Extensive  circulation  and  mixing  is  known  to  take 

place  with  n  clouds,  in  addition  to  over-all  movement  and  structural  rearrange- 
(24)  '251 

ments.  '  '  Both  directional  and  turbulent  currents  of  very  dense  water 

vapor  suggest  a  "tighter"  correlation  than  that  found  in  free  air.  Records  of 
refractometer  flights^^^^^^^^  through  clouds  substantiate  this,  and  indicate  that 
the  scale  of  turbulence  is  probably  of  the  order  10  to  20  feet.  These  same 
records  suggest  that  RMS  excursions  of  the  index  of  refraction  may  be  as  large 
as  10  or  20  N  units.  Such  results  characterize  a  turbulent  medium  far  more 
active  and  influential  for  microwave  frequencies  than  was  considered  in  Part 
II.  We  shall  therefore  undertake  a  separate  determination  of  propagation 
through  cloud  structures. 


727) - 

'  'J.  S.  Malkus  and  C.  Ronne,  Wood't  Hole  Oceanographic  Institution,  Ref 

No.  54-18,  March,  1954. 

S.  Malkus  and  R.  S.  Scorer,  Wood's  Hole  Oceanographic  Institution,  Ref. 
No  54-5,  January,  1954. 

M.  Crain,  A.  P.  Deam,  and  J.  R.  Gerhardt,  Proc.  I.R.E.,  253 

(1953). 

G.  Birnbaum  and  H.  Bussey,  Journal  of  Re  sear  ch ,  51,  171  (1953). 
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A.  Cloud  Scattering 


We  may  estimate  the  effect  of  (uniform)  turbulence  scattering  in  a 
cloud  layer^  of  thickness  L  with  the  results  developed  in  Part  II.  Consider  a 
cloud  at  mean  distance  R  from  the  receiver  as  shown  in  Figure  18.  Expressioi 
(2.  11)  yields : 
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1  •  cosp 


(4.1) 


We  shall  return  to  a  numerical  evaluation  of  this  expression  for  the  relevant 
cloud  turbulence  parameters  in  Part  C. 


B.  Direct  Wave  in  Clouds 
(1)  Phase  Variation 

We  shall  exploit  the  general  method  of  Section  II  B  to  calculate  the 
phase  shii.  n  the  direct  ray  due  to  turbulence  within  a  cloud.  It  will  be 
assumed  that  the  turbulence  pattern  within  the  cloud  is  stationary  with  respect 
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to  the  cloud.  The  line-of-aight  sweeps  through  this  pattern,  either  because  of 
the  cloud's  over  all  motion  or  due  to  an  actual  translation  of  the  propagation 
path.  The  initial  line-of-sight  is  denoted  by  and  t  seconds  later  by 
The  space  interval  between  the  two  points  and-/^  is  d  =  vt  where  v  is  the 
velocity  of  the  cloud  and/or  path  (see  Figure  19). 

The  autocorrelation  of  the  phase  shift  is  computed  as  in  Eq.  (2.20), 
with  the  correlation  (space)  interval  defined  by: 


r,2  =  +  . 


(4.2) 


Therefore; 
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The  change  of  variable  y  =  1  + - —  gives: 
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If  hi 1,  the  upper  limit  may  be  relaxed  to  infinity  with  small  error  and 
the  second  term  neglected  in  comparison  with  the  first.  The  integration  may 
then  be  performed  analytically. 


8tt^  10'^^<AN^>/ 

<a(t)a(t  +  t)>^ - 2 

.2 


(4.5) 


where  is  the  modified  Bessel  function  of  the  second  kind.  One  may  verify 
that  this  expression  reduces  to  Eq.  (2.  19)  when  t  =  0.  The  cosine  transform  of 
£q.(4.5)  yields  the  power  spectrum 


8it^lO‘*VAN^>/  /  1.  V 

W(f,  = - 5 - S(_} 


(4.6) 


which  is  plotted  in  Figure  20. 


A  ray  passing  through  a  turbulent  cloud  is  perturbed  in  two  ways  in 
that  it;  (1)  is  bent  on  entering  and  again  on  leaving,  and  (2)  suffers  a  random 
deviation  of  its  direction  of  propagation  due  to  the  cumulative  effect  of  the 
turbulence  within  the  cloud.  These  two  effects  will  cause  the  apparent  source 
seen  through  the  cloud  to  wander  about;  much  as  a  source  of  (visual)  light 
viewed  through  an  inhomogeneous  window  pane  apparently  changes  position  as 
the  line  of  sight  traverses  it. 

Figure  21  pictures  an  idealized  sharp  boundary  for  a  typical  cloud. 
The  angle  of  incidence  0.  and  refraction  0^  are  measured  with  respect  to  the 
cloud's  local  normal  to  the  boundary.  Snell's  law  relates  these  angles. 
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n.  sinO.  =  n  sin6 
1  1  o  o 


(4.7) 


Since  the  effect  of  diffuse  boundaries  is  quite  small,  one  may  put 

=  (4.8) 

where  ^  is  the  (small)  deviation  produced  at  the  boundary,  and 

n  '  =  n.  +  AN  10“^  (4.9) 

Since  AN  cr*  10.  may  expand  Eq.  (4.  7)  to  the  first  order  in  AN  and  . 

=  AN  10“^  tane.  =  AN  lO'^e.  .  (4.10) 

Here  we  have  recognized  that  n.  is  very  nearly  unity.  If  we  assume  that  the 
cloud  moves  through  the  (stationary)  line*of-sight  with  velocity  v,  the  auto¬ 
correlation  function  for  y  is 

+  t)>  =<  AN(t)  AN(t  +  t)>  +  t)>  ,  (4.11) 
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(2.  12) 


where  p(r)  is  given  by  Eq.(2.2)  and  R(r)  is  the  normalized  autocorrelation  of 
the  angle  of  incidence  0.  . 

We  know  very  little  about  the  detailed  appearance  of  cloud  boundaries 
(i.e.,  0.)  so  that  we  assume  for  convenience  that  R(r)  =  p(r).  We  shall  show 
later  that  this  choice  has  little  significance  in  considering  the  total  deflection 
caused  by  the  cloud  since  the  interior  dominates  over  the  boundaries.  The 
spectrum  for  due  to  the  fluctuating  boundary  is  thus: 
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The  deviation  produced  by  turbulence  within  the  cloud  must  be  added 
to  Eq.  (4.  13).  We  choose  the  situation  shown  in  Figure  19  and  tl  ;  rectangular 
coordinates  oriented  with  the  z*axi8  along  the  initial  line-of-sight  and  the 
x-axis  along  the  direction  of  cloud  motion.  The  ray  deviations  are. 


and 


(4.  14) 
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so  that  the  corresponding  aul  ^correlations  become 
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Using  Eq.(2.42)  and  the  rectangular  representation  for  the  interval  r^^  * 


/ - Z - Z - T 

■•12  =  y  (>‘2  -  ’‘i*  +  (yz  •  yp  *  <"^2  *  ■ 

one  can  average  the  gradients  in  Eq.  (4.16)  and  Eq.  (4.  17)  as, 


1  dpN^  1  dp  "i 

*’12  ^12 


(4.  18) 


(4.  19) 


where  u  =  x  or  y,  and  ^  =  x^  -  or  y2  -  y^  respectively.  We  shall  actually 
compute  Eq.  (4.  19)  for  ^2  '  ’^l  “  ^2  ”  ^1  “  expression  (2.2)  for 

p(r). 
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The  integration  techniquen  of  Section  11 B  may  be  applied  here  to  yield: 
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The  cosine  transform  of  this  expression  gives  the  spectrum  as: 
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which  is  plotted  in  Figure  22. 

Using  the  same  procedures,  one  obtains  from  Eq.  (4.  17) 


(4.23) 
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the  spectrum  of  which  is  also  plotted  in  Figure  22. 
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The  remarks  made  in  Part  II B  concerning  the  divergence  of  integrals  over 
these  spectra  apply  here  also. 

Note  that  the  coefficients  of  Equation8(4.  1 3),  (4.22),  and  (4.24)  all 
contain  factors  which  are  of  the  same  order  of  magnitude  except  that  Eq.  (4.  13) 
involves  while  the  others  have  the  factor  L.  We  have  consistently  assumed 
that  L  so  that  it  is  primarily  the  interior  of  the  cloud  and  not  the  refraction 

at  the  boundaries  which  contributes  to  the  deflection  of  the  ray. 


C.  Numerical  Estimates 


Cloud  parameters  are  conspicuous  by  their  absence  in  the  experimental 
literature.  To  illustrate  the  foregoing  results  however  we  shall  assume  - 
20  feet,  L  =  5000  feet,  X  =  0.2  feet,  and! 


1^AN^> 


=  10.  The  mean  square 


phase  deviation  due  to  scattering  is  obtained  from  Figure  4,  with  p  =  4  . 


<a^  >>  =  10’^  rad^  .  (4.25) 

^  SC' 

The  corresponding  direct  wave  result  from  Eq.  (i.  19)  is: 

,>  =  10’^  rad^  .  (4.26) 

^  direct- 

Clouds  thus  constitute  an  important  source  of  phase  deviation  for  any  propagation 
process. 
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Half  power  bandwidths  for  the  spectra  can  be  calculated  as  before. 

Two  different  physical  situations  present  themselves.  One  can  imagine  the 
line-of-sight  moving  through  a  cloud  because  of  the  tracking  of  a  fast  moving 
target.  In  this  case  a  velocity  of  500  fps  is  quite  reasonable.  With  the  above 
assumptions  this  yields  a  bandwidth  of  about  3  cps.  On  the  other  hand,  if  we 
consider  a  stationary  line-of-sight  with  a  cloud  drifting  through  at  v  =  20  fps. 
one  obtains  a  bandwidth  of  about  0.  12  cps.  The  preceding  estimates  apply  to 
phase  deviations  only.  Bandwidths  for  the  angular  deviation  will  be  comparable. 
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